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ABSTRACT 

We present an analysis of a 50 ks XMM-Newton observation of the merging galaxy cluster 
C1J0152.7— 1357 at z — 0.83. In addition to the two main subclusters and an infalling group detected 
in an earlier Chandra observation of the system, XMM-Newton detects another group of galaxies 
possibly associated with the cluster. This group may be connected to the northern subcluster by a 
filament of cool (1.4lg ^ keV) X-ray emitting gas, and lies outside the estimated virial radius of the 
northern subcluster. The X-ray morphology agrees well with the projected galaxy distribution in new 
K-band imaging data presented herein. We use detailed spectral and imaging analysis of the X-ray 
data to probe the dynamics of the system and find evidence that another subcluster or group has 
recently passed through the northern subcluster. C1J0152.7— 1357 is an extremely dynamically active 
system with mergers at different stages occurring along two perpendicular merger axes. 
Subject headings: cosmology: observations - galaxies: clusters: general - galaxies: high-redshift galax- 
ies: clusters: individual: (C1J0152.7— 1357) - intergalactic medium - X-rays: galax- 
ies 



1. INTRODUCTION 

Clusters of galaxies are believed to form hierarchically 
via the collapse and merger of smaller structures. Mas- 
sive clusters form most recently in this scenario, and in- 
deed, X-ray substructure which is indicative of merger 
and formation activity is f ound to be more pre valent in 
clusters at high redshifts Ipeltema et alJ 120051) . Obser- 
vations of merging systems at high redshift thus enable 
the study of the process of cluster formation. X-ray ob- 
servations provide a useful tool in this study, enabling 
measurements of the properties of the dominant bary- 
onic component of clust ers and detecting merg er-related 
features lik e shocks ('e.(/. |Markevitch et all 2002D and cold 
fronts f^.lMarkevitch et a,Ul2000llVikhlinin et al 12001(1 
in the X-ray emitting gas. 

In cosmological simulations, the most massive clusters 
form at the intersection s of the filaments of structure 
(e.g. I.Ienkins et aTlll998D . While filamentary structures 
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are observed in the larg e scale distributions of galaxies 
{e.g. iCclless e^ alJ|2Q(Jl f). thev are rarely observed in X- 
rays jScharf et al.l2000tlDurret et alJ2003HEbeling et all 
12004^ because of their low gas densities. Cluster merg- 
ers have also been detected in X-rays at the meet- 
ing points of apparent filaments, for exa mple in Abell 
85 iDurret et all 12001 [200l . Abell 52 1 llArnaud et alj 
2QQft lFerr^ttT^t*ajT2005 ) and Coma (e.g. \Ne\unsmxT^li\\ 
2003). C1J0152.7— 1357 is an example of a massive, high- 
redshift merging system located in a network of large- 
scale structures l|Kodama et al.l l2005L and so presents a 
valuable opportunity to observe the hierarchical forma- 
tion of structure. 

C1J0152.7— 1357 was discovered i n the Wide Angle 
ROSAT Pointed Survey (WARPS: IScnarf et a.l J IT997I: 
lEbeling et all 1200(1 . It was also discov ered indepen- 
dently in the RDCS lIRosati et alJ Il99l and SHARC 
ijRomer et ahl fcOOO) surveys and was sp ectroscopically 
confir m ed at a redshift of z = 0.83 lEbeling et alJ 
(2000); iDella Ceca eTall pOOOj) . Analysis of the dis- 
covery ROSAT data found the system to be highly lu- 
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minous with considerable X- ray substructure con sisting 
of two probable subclusters l)Ebeling et al J 12000). The 
cluster has also been the subj ect of BeppoSAX obser- 
vations iDella Ceca et al] |2000l) and Sunyaev-Zel'dovich 
effect imaging Ipov et alJl200l . A Chandra ACIS-I ob- 
servation enabled the two subclusters to be fully resolved 
and detected a galaxy group near the cluster to the east 
ijMauehan et al.ll2003|) . The Chandra data were used to 
measure the temperatures of the two subclusters sepa- 
rately for the first time, confirming that they are both 
hot (~ 5.5 keV) and massive. 

More recently, C1J0152.7— 1357 has bee n the subject of 
detaile d studies at optical wavelengths. iDemarco et alJ 
( 2005) confirmed 102 cluster member galaxies with VLT 
spectroscopy and found that the eastern group was at 
the same redshift as the cluster. iGirardi et al. (2005) 
used this redshift information to perform a detailed dy- 
namical analysis of the system. The distribution of 
galaxies around this system was probed out to scales 
of ~ 5 Mp c using photometric r edshifts obtained with 
Subaru bv iKodama et all {2005). In addition, this re- 
markable system has been t he subject of we ak lensing 
analyses, both ground ba sed l)Huo et alJl2004fl and with 
the HST (|Jee et al.ll2005|) . Broadly speaking, the distri- 
butions of X-ray emitting gas, cluster galaxies, and dark 
matter are found to be similar. 

In this paper we present the results of an analy- 
sis of an XMM-Newton observation of C1J0152.7-1357 
along with new K-band imaging. The large collect- 
ing area of XMM-Newton enables the most detailed X- 
ray study yet of this system. A ACDM cosmology of 
H = 70 km s- 1 Mpc -1 , and fl M = 0.3 (fi A = 0.7) 
is adopted throughout, and all errors are quoted at the 
68% level. At the cluster redshift of 0.83, 1" corresponds 
to 7.6 kpc in our chosen cosmology. 

2. DATA REDUCTION 

2.1. XMM-Newton data 

C1J0152.7-1357 was observed by XMM-Newton for 
50 ks on 2002 December 24 (ObsID 0109540101). The 
data were reduced and analysed using the XMM-Newton 
Science Analysis Software (SAS) version 6.1 with the lat- 
est calibration products available in April 2005. A binned 
image of the field observed by XMM-Newton EPIC is 
shown in Fig. ^ The target cluster and the galaxy NGC 
0720 (detected by chance at the edge of the field) are 
labelled. 

Periods of background flaring were detected and re- 
moved by applying an iterative sigma-clipping algorithm 
to lightcurves of the data from the PN and two MOS de- 
tectors. This was first done for lightcurves of the whole 
field in the 10 — 15 keV band with time bins of 50 s and a 
3cr threshold. The observation was virtually unaffected 
by background flares, and after cleaning 48 ks (53 ks) 
of useful PN (MOS) data remained. The results of this 
cleaning were checked by applying the same procedure 
to lightcurves of source- free regions in the 2 — 15 keV 
band with 200 s bins and a 2a threshold. The results of 
this stricter broad-band cleaning were completely con- 
sistent with those of the high-energy band cleaning, and 
the good time intervals defined in the 10 — 15 keV band 
were applied to the data used in all further analysis. 

In the spectral analysis of extended sources, the precise 
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Fig. 1. — Binned X-ray image of the field of view. Labels indicate 
the target cluster and the galaxy NGC 0720 which was observed 
by chance at the edge of the field. The annulus indicates the region 
used for the measurement of the background emission. 



modeling or subtraction of the background emission is 
crucial, particularly in regions of low surface brightness. 
The compact angular size of C1J0152.7— 1357 relative to 
the XMM-Newton field of view means that the back- 
ground emission local to the source can be used obtain 
background spectra. The annulus used for measuring the 
local background emission is marked in Fig. ^ and has 
inner and outer radii of 7' (3.2 Mpc) and 10' (4.6 Mpc) 
respectively. Note that all sources within this region were 
excluded in our analysis, but these are not marked on the 
Figure in order to preserve clarity. As this background 
region is further from the optical axis than the source, 
the emission in the region will be vignetted (the mean 
effective area in the background region is ~ 60% of that 
at the optical axis). The effects of vignetting were ac- 
counted for in the imaging analysis with the use of expo- 
sure maps, and in all spectral analysis by using the SAS 
task evigweight along with on-axis effective area files. 

All spectral analysi s was also perfo r med using the 
blank sky data sets of iRead fc Ponm anl i200 3l) with the 
"double subtraction" method o f lArnaiid et aTTll27)02T) . In 
this method background spectra are extracted from the 
same detector regions as the source spectra, but from the 
blank sky datasets. A background region of the source 
dataset (in this case the annulus marked in Fig. Q is 
used to derive a residual spectrum to account for any 
differences (predominantly at low energies due to dif- 
ferent levels of soft Galactic background) between the 
source and blank sky fields. The final background spec- 
trum consists of the blank sky spectrum plus the residual 
spectrum (scaled for any differences in extraction area). 
The exposure times of the blank sky observations were 
normalised so that the count rates of events detected out- 
side the telescopes' field of view in the target and blank 
sky datasets matched. Fig. |21 shows the PN background 
spectra extracted from the background annulus region in 
the source and blank sky datasets. The spectra agree well 
in the fluorescent lines at 1.5 keV and 8 keV but the nor- 
malisation of the continuum in the blank sky spectrum 
is too high. This indicates that the level of the particle- 
induced background (which is responsible for the out of 
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Fig. 2. — XMM-Newton PN background spectra obtained from 
the background annulus region in the source field (black) and blank 
sky dataset (grey). The blank sky spectrum was normalised by the 
flux of events detected outside the field of view. 



field of view counts and dominates the fluorescent line 
flux) is higher in the source dataset than the blank sky 
observations. Alternatively if the blank sky spectra were 
normalised to match the continuum level, the normalisa- 
tion of the fluorescent lines would be incorrect. Both of 
these differences would be compensated for to some ex- 
tent by the residual spectrum that is added to the blank 
sky spectrum. 

Despite these differences, almost all of the tempera- 
tures measured with the double subtraction blank sky 
background method were consistent with those measured 
by adopting a local background alone (see M5.3I for the 
exception). This is most likely due to the correcting 
effect of the residual spectrum, and the fact that the 
slope of the blank sky background spectrum continuum 
is the same as the local background. It should be noted, 
however, that the statistical uncertainties on the spec- 
tral properties measured in the low surface brightness 
regions of this data where the background contribution 
is important are large. Thus the effects of the systematic 
differences between the local and blank sky backgrounds 
cannot be sensitively tested. With those caveats, we con- 
clude that a simple local background, extracted from the 
source dataset, and vignetting-corrected with evigweight 
is the more reliable method for this particular dataset, 
and this is used for all spectral results presented here. 

2.2. Near-infrared data 

K band observations of C1J0152. 7-1357 were made 
with IRIS2 (the near-infrared imager and spectrograph) 
on the 3.9m Anglo-Australian Telescope, Siding Spring, 
Australia, on 2003 September 5, and 2004 January 10. 
The conditions on both occasions were clear, and the 
seeing was w 1.7" and 1.3 in September and January 
respectively. 

The images were dark-subtracted, flat-fielded and mo- 
saiced using standard techniques with the iraf software 
package. Flatfields were created by median combining 
the jittered object frames. The total exposure time of all 
observations was 205 minutes. 

3. IMAGING ANALYSIS 



A spectrally-weighted exposure map was produced for 
each EPIC camera, to account for the energy dependence 
of the telescope vignetting function. Exposure maps were 
produced in narrow energy bands, within which the vi- 
gnetting function varied little, and these were weighted 
according to the relative contribution of each band to an 
assumed spectral model. An absorbed MeKaL model 
with kT = 5.5 keV appropriate for C1J0152-7-1357 
(|Maughan et al.ll2003l and gj), was used. All of the X- 
ray imaging analysis was performed in the 0.3 — 5 keV 
energy band. 

A mosaiced image of the emission detected by the three 
cameras was then produced, and exposure corrected. 
This image was adaptively smoothed to show real fea- 
tu res detected at the 3c r level with the asmooth algorithm 
of Ebcling et alJ lj2005D . and logarithmically spaced con- 
tours of this smoothed emission are shown overlaid on a 
NIR image in Fig. [3] All of the sources in Fig. |3| except 
for those labelled are point sources, with the possible ex- 
ception of the north eastern most source at a[2000.0] = 
01 h 52 m 45.64 s <5[2000.0] = -13°55'29.0" Base d on 
the redshifts measured by iDemarco eta!] ©05), the 
brightest X-ray source near the subclusters (a[2000.0] = 
01 h 52 m 43.62 s <5[2000.0] = -13°59'1.0") appears associ- 
ated with a cluster member galaxy, and the point source 
at a[2000.0] = 01 h 52 m 34.62 s <5[2000.0] = -13°59'30.2". 
Within the size of the XMM-Newton point spread func- 
tion (PSF), there are two cluster member galaxies and 
one foreground galaxy coincident with the point source 
at a[2000.0] = 01 h 52 m 39.75 s 6[2000.0] = -13°57'42.8". 

The complex morphology of the extended emission in 
C1J0152.7— 1357 is immediately apparent in Fig. [3J The 
north and south subclusters (labelled C and E in Fig. 
EJ) and the group to the east (D) were all detected in 
the C handra observation of this system ijMaughan et alJ 
2003). In addition, XMM-Newton also detects an en- 
velope of low surface brightness emission surrounding 
the subclusters and eastern group, which extends in a 
apparent filament to a second group (A) to the north- 
west. The redshift of this group is unknown. However, 
while the associated galaxies were not inc luded in the 
VLT observations of D emarco et aj] {2005), the group s 
position appears to coincide with a clump of galaxies 
of the same photometr ic redshifts as C1J0152.7— 1357 
in iKodama et alJ l)2005|) . We address the question of 
whether this north-west group is associated with the 
main cluster in more detail later. The morphology of 
the system suggests intersecting north-east to south-west 
and south-east to north-west merger axes, which are in- 
dicated in Fig. |3J 

Fig. also shows that the X-ray emission from the 
north subcluster is extended in the northern direction, 
with a possible second X-ray peak (labelled B). This ex- 
tension is in the direction of an overdensity of galaxies 
slightly further north, which als o appears as a peak in 
the weak lensing mass map of Dee et alJ l|2005|) . 

For convenience, the labelled regions in Fig. [3]are sum- 
marised in Tabled 

The distribution of NIR light in C1J0152.7-1357 was 
also investigated with the K band imaging data. The 
NIR luminosit y of a galaxy is well c orrelated with its dy- 
namical mass l)Gavazzi et al.lH"99(^ . making the K band 
an excellent choice with which to trace out the distri- 
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Fig. 3. — Adaptively smoothed contours of the X-ray emission detected from C1J0152-7— 1357 with XMM-Newton are overlaid on the 
K band IRIS2 image. The dashed lines indicate the two apparent merger axes of the system and the labels are described in the text and 
summarised in Table HI 



TABLE 1 

Summary of the structures labelled in Fig. |3] 



Region 


a[2000.0] 


<5[2000.0] 


Description 


A 


01 h 52 m 34.18 s 


-13°54'42.3" 


North west group of unknown redshift 


B 


01 h 52 m 45.15 s 


-13° 57' 02.4" 


Northern extension of north subcluster 


C 


01 h 52 m 43.76 s 


-13°57'24.4" 


North subcluster 


D 


01 h 52 m 51.99 s 


-13°58'13.2" 


Galaxy group at cluster redshift 


E 


01 h 52 m 39.83 s 


-13°58'26.7" 


Southern subcluster 



bution of galaxies within the cluster. The image was 
divided into a grid of 1000x1000 bins, and the number 
of galaxies detected at a significance > 3cr in each bin was 
counted. Objects which were brighter than the brightest 
known members were excluded. Due to the poor seeing, 
the galaxy density may be underestimated in the high- 
est d ensity regions. Howeve r, all of the known members 
from iDemarco etall (2005) were included by manually 
dcblcnding some of the galaxies. An image of this map 
was created and adaptively smoothed at the 90% level. 
Contours of galaxy number density are shown overlaid 
on a K band image and a smoothed X-ray image in Fig. 



0J There is a clear similarity in the galaxy and X-ray 
distribution; both the subclusters and groups are appar- 
ent. An overdensity of galaxies to the north of the X-ray 
emission in region B is also clear. 

The comparison of Figs. |3] and 0] shows a clear off- 
set between the peak of the galaxy light and X-ray 
emission in the southern subcluster. This was first no- 
ticed in the Chandra observation of C1J0152.7— 1357 by 
iMauehan et alJ l)2003j) who suggested that the offset was 
due to the collisionless galaxies (and presumably dark 
matter) moving ahead of the X-ray ga s which is slowed 
by ram pressure due to the merger. Dee et all (2005) 
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found similar offsets between the mass peaks of both 
subclusters in their weak lensing analysis and the X-ray 
centroids, which is consistent with this explanation. We 
note, however that the X-ray peak of the northern sub- 
cluster is coincident with that of the galaxy distribution, 
but the northern extension of that subcluster (likely due 
to another merger) leads to an apparent offset between 
the X-ray centroid and galaxy distribution. 

3.1. X-ray surface brightness modeling. 

The two-dimensional (2D) X-ray surface brightness 
distribution in C1J0152-7— 1357 was modeled in Sherpa 
with an elliptical (3— model for each of the two subclusters 
and two groups, and a fla t and a vignetted background 
(see iMauehan et alJ 12004 for a more detailed descrip- 
tion of the method). Due to the low surface brightness 
of the two groups, their model ellipticities were fixed at 
zero and their core radii were fixed at 100 kpc. This 
core radius is appropriate for the temperatures o f these 
groups as measured in $11 (San derson et alJ l2003). The 
models were fit simultaneously to binned images for each 
EPIC camera after being convolved with the appropriate 
PSF and exposure map. All model parameters except for 
the amplitudes were tied between the different detectors. 
The background levels were determined from a fit to a 
local source free region and fixed. 

The data residuals from the best-fitting models for each 
camera were combined and smoothed with a Gaussian 
of a = 12" (3 image pixels). The resulting image is 
shown in Fig. Point sources were excluded both dur- 
ing the fitting process and from these residuals. There 
are several regions of excess emission above the model 
visible in Fig. We note that some of these resid- 
ual features occur near or across PN CCD gaps. How- 
ever, they are also present in the residuals for the MOS 
detectors alone so are not due to exposure correction 
problems around the CCD gaps. In Fig. |HJ contours 
of the positive surface brightness residuals are overlaid 
on the best-fitting surface brightness model, and on a 
NIR image. In this Figure, the surface brightness mod- 
els for each detector have been summed, after convolu- 
tion with the PSF and multiplication by the exposure 
map. The PN CCD gaps are thus visible in the surface 
brightness model. The same labels used in Fig. arc 
included to aid orientation. Positive residuals are asso- 
ciated with the subcluster cores (C and E), which may 
indicate the presence of dense cool cores, and also with 
region B. The region of excess emission at a[2000.0] = 
01 h 52 m 40.75 s (5[2000.0] = -13°56'33.1" to the west of re- 
gion B is aligned with several faint galaxies, while that 
at a[2000.0] = 01 h 52 m 42.71 s 5[2000.0] = -13°59'52.5" to 
the south east of the southern subcluster (E) coincides 
with several brighter galaxies. Both features appear to be 
associated with galaxies whose photometric redshifts are 
consistent with that of C1J0152.7— 1357 as measured by 
iKodama et alJ l)2005|) . and both appear in the contours 
of galaxy light in Fig. 0] While this is not statistically 
rigorous, it is likely that these regions of excess emis- 
sion are linked to levels of substructure in the system 
which are not present in our simple surface brightness 
model. Similarly, the negative residuals to either side of 
the northern subcluster core in Fig. 0are probably due 
to the inadequacy of a simple elliptical model to describe 
the subcluster. 



There is excess emission between the the two main 
subclusters (C and E). This excess was also detected 
in the Chandra obs ervation of C1J0152.7— 1357 by 
IMauehan et alJ (12003ft . A possible explan ation for this 
feature proposed bv lMaugharTet al. ( 2003) is a region of 
increased density due to the compression of the gas in the 
merger. Howeve r, the weak lensi ng mass reconstruction 
of this system bv lJee et all (|2005|) shows two mass peaks 
in this region, and there are also several cluster galax- 
ies in this region (see Fig. 0). This suggests that the 
excess emission may be due to gas associated with mass 
clumps in this region, which may be additional subgroups 
or clusters in this complex system. 

3.2. Projections of the X-ray surface brightness 

An alternative way of visualising the data is to use 
projections of the surface brightness along different axes. 
In the case of C1J0152.7— 1357 the two merger axes are 
natural choices. A rectangular region aligned with each 
of the merger axis was defined. The X-ray counts in a 
combined PN and MOS 0.3 — 5 keV image within this 
region were then projected onto its long axis. The same 
procedure was applied to combined PN and MOS image 
of the best-fit 2D model, after convolution with the PSF 
and exposure map. The data projection was adaptively 
binned to give fractional uncertainties of < 10% on each 
point. 

Fig. [3 shows the projection along the south-east to 
north-west merger axis. For this axis, a region of width 
75" centred on the two groups was used. Several fea- 
tures are labelled in Fig. using the same labels as 
Fig. |3and Tabled The region between C and A cor- 
responds to the apparent filament in Fig. |3] In this 
region, the model represents the simple superposition of 
the group and cluster emission and the data show no 
significant departure from this model. The filamentary 
appearance of the emission in this region could thus sim- 
ply be due to the superposition of the cluster and group 
emission. However, given the difficulty of modeling the 
surface brightness distribution of this system, we cannot 
rule out the possibility of filamentary emission in this 
region. The point source visible in Fig. halfway be- 
tween B and A was excluded from the data and model 
in this projection. The projection of the north-east to 
south-west merger axis is shown in Fig. |H1 A rectan- 
gular projection region of width 50" centred on the two 
subclusters was used. The substructure around region B 
in Fig. |31 is visible in this projection. Excess emission be- 
tween the subclusters (C and E) is apparent, and a bright 
point source which was excluded from the 2D modeling, 
but not from the data in this projection is visible at the 
right end of Fig. |H1 The disagreement between the model 
and data at the position of the chip gap at 230" is due 
to the imperfect exposure correction around the chip gap 
combined with the bright point source. 

4. SPECTRAL MAPPING 

The high quality X-ray data enable spatially resolved 
temperature measurements of the gas in this merger. 
A met hod based on that described bv 10 'Sullivan et all 
( 2005) was used to produce the temperature map. First, 
a "radius map" was produced, based on a binned X- 
ray image. The radius map recorded the radius of the 
circular region enclosing > 1000 background-subtracted, 
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Fig. 4. — Left: Contours of constant galaxy number density derived from the AAT K-band images are overlaid on the same NIR image. 
The contours have been adaptively smoothed at 90% significance. Right: The same contours are overlaid on the smoothed X-ray image 
from which the X-ray contours in Fig. |2]were taken. 




Fig. 5. — Smoothed X-ray image of the residuals from the best- 
fitting model to the surface brightness distribution in C1J0152- 
7—1357. Smoothed contours of X-ray surface brightness are over- 
laid. 



point-source-excluded photons centered on each pixel 
in the binned image. A maximum radius of 60" was 
imposed, and pixels which failed to meet these crite- 
ria were excluded from the radius map. Spectra were 
then extracted for each pixel within the appropriate ra- 
dius, and the three EPIC spectra were fit simultane- 
ously with an absorbed MeKaL model. During the fit- 
ting process, the model parameters were tied except for 
the normalisations which were independent. The metal 
abund ance was fixed at 0.3Z(n u sing the abundance ta- 
bles of lAnders fc Grevessd l)1989l) , an d the absorbing col- 
umn was fixed at the Galactic value l)Dickev fc Lockmanl 
1990). All spectra were grouped to have at least 20 total 
counts per spectral bin, and were fit in the 0.3 — 7 keV 
band using local background spectra extracted from the 
background annulus region (Fig. [Q. 

The resulting projected temperature map is shown in 
Fig. E| It should be noted that the method used for pro- 
ducing the temperature map means that adjacent pixels 



are not independent. The radius map in Fig. E| gives 
the extent of the spectral extraction region used at each 
pixel. Broadly speaking, the gas in the system has tem- 
peratures in the range 4 — 7 keV but the spatial distri- 
bution of temperatures is not smooth. There are three 
apparently hotter regions in the temperature structure 
which are labelled in Fig. El The northern-most (kTl) is 
aligned with the residual emission at region B (see Figs. [S] 
andlHI) • The hotter regions at kT2 and kT3 are not associ- 
ated with any surface brightness features. The mixing of 
spectra between adjacent pixels in this spectral mapping 
method can make such features appear more significant 
than they truly are, as they occur in several adjacent pix- 
els. The temperatures and fit parameters of the hottest 
pixel in each of these regions are given in Tabled While 
region kTl is a local maximum in temperature, it is ac- 
tually no hotter than the average temperature across the 
cluster of 5.5 keV. Regions kT2 and kT3, however, are 
hotter at the ~ 2a level. With the exception of pixels 
close to the boundaries of the temperature map, none of 
the measured temperatures are significantly cooler than 
5.5 keV. 

A simple alternative method of producing a temper- 
ature map was also implemented. In this method geo- 
metrical regions were defined by hand using the surface 
brightness distribution and residuals as a guide. Spectra 
were extracted from these regions (with point sources 
excluded) and fit as above. Again, vignetting-corrected 
background spectra extracted from the background an- 
nulus region were used. The goal in following this method 
was to isolate some of the features in the system and ob- 
tain spectra that are free from the mixing effects of the 
previous method. We were also able to extract and fit 
spectra for regions which did not pass the binning algo- 
rithm's criteria. Fig. 1 101 shows the resulting temperature 
map and the fit parameters of the numbered regions are 
given in Table El The table also gives the net source 
counts detected with each EPIC camera. We note that 
three of the regions (1, 2 and 7) have low net counts 
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Right ascension 

Fig. 6. — Left: The best-fitting two-dimensional surface brightness model of C1J0152.7— 1357 with contours of the positive residuals 
overlaid. Right: The same contours are overlaid on an IRIS2 K band image, with the same labels used in FiglSI 



data 




50 100 150 200 250 300 350 400 450 
Distance (arcsec) 

Fig. 7. — Projection of the surface brightness in a 75" wide 
region onto the south-east to north-west merger axis. The points 
show the projection of the imaging data and the solid line is the 
projection of the best-fit 2D model. The dashed line shows the 
locally estimated background level. 




Distance (arcsec) 

Fig. 8. — Projection of the surface brightness in a 50" wide 
region onto the north-east to south-west axis. The points show the 
projection of the imaging data and the solid line is the projection of 
the best-fit 2D model. The dashed line shows the locally estimated 
background level. 



TABLE 2 

Values of the spectral fit parameter s fo r the regions 
shown in Fig. |3] and Fig. fTTJl 



Net Counts 

Region kT (keV) ~TN MOS1 MOS2 x 2 » 



kTl 


^ n+0.9 
°- y -0.6 


579 


240 


276 


50.8 


59 


kT2 


10.3+?° 


549 


244 


245 


46.0 


52 


kT3 


7 o+1.4 


705 


250 


228 


44.1 


59 


1 


1 5 +U ^ 


89 


31 


38 


18.5 


12 


2 


i-4±S:? 


150 


18 


28 


102.8 


63 


3 




340 


100 


123 


48.4 


51 


4 


5 6+ 11 
°- D -0.8 


355 


108 


135 


20.6 


27 


5 


6 4+°- 8 


625 


275 


233 


52.0 


53 


6 


5 4+°- 4 
°- 4 -0.4 


1425 


620 


682 


157.7 


159 


7 


2 4+°- 6 


179 


76 


61 


26.5 


31 


8 


t- n+0.8 
J.O_q 7 


506 


161 


170 


42.6 


44 


9 


4 Q+0-6 


499 


188 


212 


45.4 


38 


10 


e. Q+0.4 
O.o_ 4 


1590 


492 


565 


142 


144 



(150 — 300). This means that the grouping of the spectra 
to 20 total counts per bin and the use of the % 2 statistic 
may not be appropriate. In these cases, ungrouped ver- 
sions of the spectra were also fit using the C statistic in 
XSPEC. In all cases, the temperatures and uncertainties 
were in good agreement with the ones quoted in Tabled 
This method indicates that the two main subclusters 
and their surrounding regions have temperatures con- 
sistent with 5.5 keV. The eastern group (region 7) is 
significantly cooler at 2.41q4 keV, with a temperature, 
appropriate for a group or small cluster. The redshift 
of the north-west group (region 1) cannot be determined 
from the X-ray spectra, but under the assumption that 
it is at the cluster redshift it too has a low temperature 
consistent with a galaxy group. The emission in region 
2, between group A and the northern subcluster also has 
a low temperature of 1.4+q'? keV assuming a redshift of 
0.83. 
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Fig. 10. — Projected temperature map of C1J0152-7— 1357 with 
contours of X-ray surface brightness overlaid. Darker regions are 
cooler and lighter regions hotter. The fit parameters of the num- 
bered regions are given in Table 121 



The regions associated with positive surface brightness 
residuals in the subcluster cores show no evidence for the 
cooler gas which would be expected if the excess emis- 
sion were due to cool cores. However, the size of the re- 
gions are small (~ 30" diameter) so the blurring effect of 
the XMM-Newton PSF would make such emission more 
difficult to detect. The Chandra observation of this sys- 
tem, with its negligible PSF, also detects positive surface 
brightness residuals in the cluster cores, but the data 
quality is insufficient to extract spectra from those re- 
gions. The statistical uncertainties on the XMM-Newton 
temperatures along with the unaccounted for projection 
effects mean that the presence of cool cores in either sub- 
cluster cannot be ruled out. 

Regions 3, 4 and 8 were also chosen based on the sur- 



face brightness residuals. These do not show departures 
from the subcluster temperatures. In particular, region 
8 between the subclusters is not significantly hotter. 

Finally, the global luminosity of the cluster was mea- 
sured. Due to the complex morphology of the system, 
extrapolating luminosities measured from spectral fits to 
the individual components of the system out to a large 
radius may lead to inaccuracies. Instead, global spec- 
tra were extracted for each EPIC camera from within a 
radius o f 1.4 Mpc (184", chosen to match the apertures 
used in iMarkevitchl Il998|) , after the exclusion of point 
sources. The resulting spectra were fit as before using 
the local background spectra, and the best fitting tem- 
perature was 3.7±0.2 keV (jfjv = 588/558). While the 
emission within this large aperture is clearly not isother- 
mal, the spectra did not support the modelling of a sec- 
ond thermal component. The bolometric luminosity (rest 
frame 0.01 — 100 keV) within this aperture was measured 
to be 1.3 ± 0.2 x 10 45 erg s _1 . The uncertainties on the 
luminosity incorporate those on the best-fit temperature 
and normalisation of the spectral model. 

5. DISCUSSION 

This deep XMM-Newton observation confirms the 
global temperatur es measured by Chan dra for the two 
main subclusters (Ma ughan et al J 120 03). and supports 
the conclusion that C1 J0152.7— 1357 is a massive merger. 
Maug han et alJ 1)20031) included just the two main sub- 
clusters in their mass analysis, but the XMM-Newton 
temperatures of the two groups can be used to esti- 
mate their contribution to the total mass of the sys- 
te m. Using the high-red shift mass-temperature relation 
of iMaughan et all l)2005|) (based on a sample of relaxed 
and unrelaxed clusters including the two subclusters of 
C1J0152.7— 1357) we estimate that an object of ~ 2 keV 
will have a total virial mass of ~ 1 x 1O 14 M . In the same 
way, the approximate temperature of the two subclusters 
(5.5 keV) gives a mass of ~ 5 x 1O 14 M0 for each. We thus 
concluded t hat the mass of the sy stem is ~ 1 x 10 15 M Q 
(in line with IMaughan et alJl2003j) and that this is dom- 
inated by the two main subclusters. 
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In the following sections we discuss the important fea- 
tures of C1J0152.7-1357 detected by XMM-Newton. 

5.1. Substructure in the northern subcluster 

The combination of the galaxy overdensity, weak- 
lcnsing mass peak and the X-ray substructure in the 
north of the northern subcluster (region B in Fig. |2J 
suggest the presence of a third subcluster or group in the 
NE-SW merger axis of the C1J0152. 7-1357. One plausi- 
ble explanation for the configuration of the galaxy, gas, 
and dark matter distribution in this region is that the 
third subcluster has recently passed through the north- 
ern subcluster, traveling in a northerly direction. The 
dark matter and galaxies of the third subcluster were 
unaffected by the encounter due to their small (or non- 
existent) collisional cross sections. Its gas, meanwhile, 
was stripped from the potential, with the northern ex- 
tension of the X-ray emission and possible second X-ray 
peak in that region being due to the surviving dense gas 
core of the third subcluster. 

The mean redshift of the galaxies associated with this 
possible third subcl uster is 0.838 ± 0.006 (using the red- 
shifts measured bv iDemarco et all 12005): the same as 
that of the northern subcluster itself within the uncer- 
tainties. This suggests that any motion is in the plane of 
the sky. 

5.2. The properties of group D 

In their dynam ical analysis of C1J0152.7— 1357, 
iGirardi et al. ( 2005) find a velocity dispersion of a v « 
700 km s . The X-ray temperature of the group 
(2.4^0 4 keV ) is slightly low co mpared to the velocity dis- 
persion (e.g. iXue fc Wull2000j) . but it is consistent with 
the a v —T relation within the scatter. Using the MeKaL 
normalisation of the spectral fit in this region, and as- 
suming that the group is spherical, the gas mass within 
a radius of 20" was estimated to be 6.0 ± 0.6 x 1O 11 M . 
This can be compared with the total mass within the 
same radius from the weak lensing data l).Iee et alJl2005l) 
to give a gas mass fraction of / gas « 0.03. At z — 0.83 the 
radius of 20" corresponds to 150 kpc, so this gas fraction 
is measured only in the central region of group D. The 
gas fraction found is consistent with that measured in the 
central regions of local systems of a similar temperature 
l)Sanderson et al.ll2003|) . 

5.3. Group A and the possible filament 

Another interesting feature of the XMM-Newton data 
is the possibility that galaxy group A is falling into the 
cluster along a merger axis which runs through the north- 
ern subcluster (C) and group D. The K-band light distri- 
bution shows similar structure to the X-ray morphology 
in this region and the galaxies associated with the group 
are faint. However it is not known if the group is at the 
same redshift as the cluster. 

The gas density in the region of emission between 
group A and the northern subcluster was estimated us- 
ing the MeKaL normalisation from the spectral fit to 
region 2 (Fig. I10|) . In order to simplify the computation 
of the volume of region 2, the MeKaL normalisation was 
first scaled up to account for the area excluded in region 
1. The volume of region 2 was then calculated assum- 
ing it is a cylinder with rotational symmetry about its 



long axis. With the simplifying assumption that the den- 
sity is constant throughout that volume, the Hydrogen 
number density in this region is 4.0 ± 0.6 x 10~ 4 cm~ 3 . 
These results should be treated with some caution as the 
emission is very low surface brightness, and the issue of 
background subtraction is critical. In particular, while 
the results above were based on spectral fits with a local 
background, we were unable to obtain acceptable spec- 
tral fits with a blank sky background. However, the un- 
certainties in accounting for the differences in soft Galac- 
tic X-ray emission in the source and blank sky datasets 
are particularly important in the study of a region of low 
surface brightness, cool gas such as this. As discussed in 
the use of a local background should to be the most 
reliable method in this instance. 

It is interesting to compare this possible filament with 
the ap parently isolated filament detected bv lScharf et alJ 
(2000) with ROSAT. The surface brightness in region 2 
is 1.3 ± 0.4 x 10~ 15 erg s _1 cm~ 2 a rcmin" 2 . This is an 
order of magnitude brighter than the lScharf et alJ (2000) 
filament. While that filament has no apparent connec- 
tion to any massive clusters, it is likely that the gas den- 
sities would be higher in regions of filaments that are 
feeding onto a massive cluster of galaxies like that in 
C1J0152.7— 1357. However, as demonstrated in i)3.2l the 
apparent filament may simply be due to the superposi- 
tion of the emission from the group and cluster. 

In order to estimate of the extent of the northern sub- 
cluster, a mass profile was derived assuming that the gas 
follows an isothermal /3-profile in hydrostatic equilibrium 
with the cluster potential. While these are clearly gross 
simplifications for this complex system, they are suffi- 
cient for the purposes of estimating an approximate virial 
radius for the northern subcluster. The gas temperature 
of region 6 in Fig. 1101 along with the /^-profile parame- 
ters from the 2D surface brightness fitting (r c = 40 ± 3", 
f3 = 0.86±0.06) were then used to estimate a virial radius 
of i?200 ~ 1-3 Mpc. Here, we define the virial radius as 
the radius enclosing a mean density that is a factor of 200 
times the critical density at z = 0.83. Numerical simu- 
lations show that this overdensity radius approximately 
separates the virialised part of clusters from the in falling 
material surrounding them l|Navarro et alJ lT995). The 
estimated virial radius of the northern subcluster falls 
approximately halfway along region 2 in Fig. Dis- 
regarding projection effects, the southern subcluster and 
group D all fall within the virial radius of the northern 
subcluster. If the group A is associated with the system, 
it is most likely at an earlier merger stage, lying outside 
the estimated virial radius. 

6. CONCLUSIONS 

C1J0152.7— 1357 is a fascinating system, comprising 
two main subclusters in an apparently early stage of 
merging and possibly two infalling groups. There also 
is evidence for late stage merger activity in the northern 
subcluster. The formation of a massive cluster is occur- 
ring along two main merger axes which appear perpen- 
dicular in the plane of the sky. This unique system pro- 
vides a dramatic example of cluster formation via merg- 
ers on different scales and at different stages, providing 
further compelling observational support for the hierar- 
chical assembly model of cluster formation. 



10 



B. J. Maughan et al. 



7. ACKNOWLEDGMENTS 

We thank Christine Jones for useful discussions of this 
work. BJM is supported by NASA through Chandra 
Postdoctoral Fellowship Award Number PF4-50034 is- 



sued by the Chandra X-ray Observatory Center, which is 
operated by the Smithsonian Astrophysical Observatory 
for and on behalf of NASA under contract NAS8-03060. 
SCE acknowledges PPARC support. 



REFERENCES 



Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Arnaud, M., Majcrowicz, S., Lumb, D., Neumann, D. M., Aghanim, 
N., Blanchard, A., Boer, M., Burke, D. J., Collins, C. A., Giard, 
M., Nevalainen, J., Nichol, R. C, Romcr, A. K., & Sadat, R. 
2002, A&A, 390, 27 

Arnaud, M., Maurogordato, S., Slezak, E., & Rho, J. 2000, A&A, 
355, 461 

Colless, M., Dalton, G., Maddox, S., Sutherland, W., Norberg, P., 
Cole, S., Bland-Hawthorn, J., Bridges, T., Cannon, R., Collins, 
C, Couch, W., Cross, N., Deeley, K., De Propris, R., Driver, 
S. P., Efstathiou, G., Ellis, R. S., Frenk, C. S., Glazebrook, K., 
Jackson, C, Lahav, O., Lewis, I., Lumsden, S., Madgwick, D., 
Peacock, J. A., Peterson, B. A., Price, I., Seaborne, M., & Taylor, 
K. 2001, MNRAS, 328, 1039 
Delia Ceca, R., Scaramella, R., Gioia, I. M., Rosati, F., & Squires, 

G. 2000, A&A, 353, 498 
Demarco, R., Rosati, P., Lidman, C, Homeier, N. L., Scannapicco, 
E., Bem'tez, N., Mainieri, V., Nonino, M., Girardi, M., Stanford, 
S. A., Tozzi, P., Borgani, S., Silk, J., Squires, G., & Broadhurst, 
T. J. 2005, A&A, 432, 381 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Durret, F., Lima Neto, G. B., & Forman, W. 2005, A&A, 432, 809 
Durret, F., Lima Neto, G. B., Forman, W., & Churazov, E. 2003, 

A&A, 403, L29, 1SS 
Ebeling, H., Barrett, E., & Donovan, D. 2004, ApJ, 609, L49 
Ebeling, H., Jones, L. R., Pcrlman, E., Scharf, C, Horner, D., 
Wegner, G., Malkan, M., Fairley, B., & Mullis, C. R. 2000, ApJ, 
534, 133 

Ebeling, H., White, D., & Rangarajan, F. 2005, MNRAS 
Ferrari, C, Arnaud, M., Ettori, S., Maurogordato, S., & Rho, J. 

2005, ArXiv Astrophysics e-prints 
Gavazzi, G., Pierini, D., & Bosclli, A. 1996, A&A, 312, 397 
Girardi, M., Demarco, R., Rosati, P., & Borgani, S. 2005, A&A, 

442, 29 

Huo, Z., Xue, S., Xu, H., Squires, G., & Rosati, P. 2004, AJ, 127, 
1263 

Jee, M. J., White, R. L., Bem'tez, N., Ford, H. C, Blakeslee, J. P., 
Rosati, P., Demarco, R., & Illingworth, G. D. 2005, ApJ, 618, 46 

Jeltema, T. E., Canizares, C. R., Bautz, M. W., & Buote, D. A. 
2005, ApJ, 624, 606 

Jenkins, A., Frcnk, C. S., Pearce, F. R., Thomas, P. A., Colberg, 
J. M., White, S. D. M., Couchman, H. M. P., Peacock, J. A., 
Efstathiou, G., & Nelson, A. H. 1998, ApJ, 499, 20 



Joy, M., LaRoquc, S., Grego, L., Carlstrom, J. E., Dawson, K., 
Ebeling, H., Holzapfel, W. L., Nagai, D., & Reese, E. D. 2001, 
ApJ, 551, LI 

Kodama, T., Tanaka, M., Tamura, T., Yahagi, H., Nagashima, M., 
Tanaka, I., Arimoto, N., Futamase, T., Iye, M., Karasawa, Y., 
Kashikawa, N., Kawasaki, W., Kitayama, T., Matsuhara, H., 
Nakata, F., Ohashi, T., Ohta, K., Okamoto, T., Okamura, S., 
Shimasaku, K., Suto, Y., Tamura, N., Umctsu, K., & Yamada, 
T. 2005, PASJ, 57, 309 

Markevitch, M. 1998, ApJ, 504, 27 

Markcvitch, M., Gonzalez, L., Vikhlinin, A., Murray, S., Forman, 

W., Jones, C, & Tucker, W. 2002, ApJ, 567, L27 
Markevitch, M., Ponman, T. J., Nulsen, P. E. J., Bautz, M. W., 

Burke, D. J., David, L. P., Davis, D., Donnelly, R. H, Forman, 

W. R., & Jones, C. 2000, ApJ, 541, 542 
Maughan, B. J., Jones, L. R., Ebeling, H., Perlman, E., Rosati, P., 

Frye, C, & Mullis, C. R. 2003, ApJ, 587, 589 
Maughan, B. J., Jones, L. R., Ebeling, H., & Scharf, C. 2004, 

MNRAS, 351, 1193 
— . 2005, MNRAS, in press, (astro 

Navarro, J. F., Frcnk, C. S., & White, S. D. M. 1995, MNRAS, 
275, 720 

Neumann, D. M., Lumb, D. H., Pratt, G. W., & Briel, U. G. 2003, 
A&A, 400, 811 

O'Sullivan, E., Vrtilek, J. M., Kempner, J. C, David, L. P., & 

Houck, J. C. 2005, MNRAS, 357, 1134 
Read, A. M. & Ponman, T. J. 2003, A&A, 409, 395 
Romer, A. K., Nichol, R. C, Holden, B. P., Ulmcr, M. P., Pildis, 

R. A., Merrelli, A. J., Adami, C, Burke, D. J., Collins, C. A., 

Metevier, A. J., Kron, R. G., & Commons, K. 2000, ApJS, 126, 

209 

Rosati, P., Delia Ccca, R., Norman, C, & Giaconni, R. 1998, ApJ, 
492, L21 

Sanderson, A. J. R., Ponman, T. J., Finogucnov, A., Lloyd-Davics, 
E. J., & Markevitch, M. 2003, MNRAS, 340, 989 

Scharf, C, Donahue, M., Voit, G. M., Rosati, P., & Postman, M. 
2000, ApJ, 528, L73 

Scharf, C, Jones, L. R., Ebeling, H., Perlman, E., Malkan, M., & 
Wegner, G. 1997, ApJ, 477, 79 

Vikhlinin, A., Markevitch, M., & Murray, S. S. 2001, ApJ, 551, 160 

Xue, Y. & Wu, X. 2000, ApJ, 538, 65 



